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Abstract 



We analyze the stochastic scaUng laws arising in the invicid limit of the decaying solutions of the Burgers 
equation. The linear scaling of the velocity structure functions is shown to reflect the domination by shocks of 
the long-time asymptotics. We exhibit new self-similar statistics of solutions describing phases with diluted 
shocks. Some speculations are included on the nature of systems whose large time behavior is described by 
the new statistics. 
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1 Introduction 



The Burgers equation, a version of the Navier-Stokes equation without pressure, takes in the 1+1 dimensions 
the form: 

dtu + udxU — vd^u = (1) 

where u — u{x, t) is the velocity field. We have not included the force term since we shall be interested 
in the free decay of initial data. Although we shall stick to the 1-dimensional space, most of the following 
can be generalized to higher dimensions. We are interested in statistical properties of the velocity field 
at time t > 0, given the statistics of random initial data. The equation has the u{t, x) i— > —u{t, —x) 
symmetry which, if respected by the statistics of the initial conditions, will persist at all times. The problem 
is to evaluate the n-point correlation functions (Wju(xj,t)) at equal time in the invicid limit v ^ Q. At 
large time these correlation functions are expected to flow towards some 'universal' functions manifesting a 
self-similar character. In other words: 

n 



t)) ~ u"it) B„(^] for t large (2) 



with l(t) and u{t) ~ dtl{t) being the characteristic length and the characteristic velocity at time t. Note 
the order of the limits: first lim and then lim . These limits do not commute. Each asymptotic universal 

u — >Q t — ^oo 

statistics, which are specified by their correlation functions Bn{xj /l{t)), have their own basin of attraction. 

If the initial statistics are Gaussian with zero mean, they are encoded in the initial velocity two-point 
function T(x — y) = (uo(a;)uo(y)) which we assume translation invariant. As is known since the work of 
Burgers, the large time behavior depends crucially on whether JJ =^ J dxT{x) vanishes or not. The case 
J was the case first studied by Burgers himself H. In that case the large time behavior is governed 
by a self-similar solution with characteristic length l{t) ^ t^^^. The case J — was analyzed by Kida in 
his important paper ||]. It is convenient to introduce the potential $(a;,i) such that u{x,t) — dx^{x,t). If 
JT" = 0, we may assume that the initial potential is Gaussian with mean zero and the two-point function 

G{x - y) = {'Poix)My)) ■ (3) 

Of course T{x) = ~dxG{x). Assuming a priori that the minima of the initial potential are independent, 
Kida showed in ref. Q that the large time behavior has a characteristic length l{t) ^ t^/^ up to a logarithmic 
correction. A more precise formulation of this statement, proved in ref. under very mild hypothesis on 
G{x), (e.g. for G{x) being a smooth function decreasing rapidly enough at infinity, c.f. also and references 
therein) states that the following limit exists: 

lim |loge|* u ^|loge|^^ UK{x,t) . (4) 

Here and in the following, = means an equality in law, i.e. inside any correlation functions. The statistics of 
the limiting velocity field UKix,t) was explicitly constructed in ||. UK{x,t) is self-similar with a diffusive 
scaling, Ixit) — i^^^ A^/"* where A = G(0). Note that one of the hypothesis for having Kida's statistics at 
large time is that < A < cx), i.e. that the initial potential two-point function is regular at the origin. 

In this letter, we shall construct other self-similar solutions of decaying Burgers turbulence. Although 
these solutions are different from Kida's statistics, they share in commun the fact to be constructed from 
Poisson point processes. These solutions may be relevant in the large time behavior of systems whose initial 
correlation functions are singular at coinciding points. 

In Sect. 2, we introduce a few basic facts concerning the Burgers equation and we describe some universal 
features of fields localized on shocks. In Sect. 3, we construct the self-similar statistics and we prove that 
they effectively are solutions of the turbulent problem. We also give a few more detailed informations on a 
particular case. Comments and speculations are gathered in Sect. 4. 

Acknov^rledgements. We would like to thank Antti Kupiainen for the collaboration in initial stages of 
this work and to Marc Mezard for stimulating discussions. 
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2 Basic facts about the Burgers equation 

In order to fix notations, we recall few elementary facts concerning the Burgers equation 0, ||. As is well 
known, the equation is solved by implementing the Cole-Hopf transformation which maps it to the heat 
equation. This works as follows. Let Z{x,t) = exp[— ^$(x, t)] where u{x,t) — dx^{x,t). Eq. for u is 
mapped into the heat equation for Z: 

[dt-vdl] Z{x,t) = Q. 

Thus, given the initial condition u{x, t — Q) = uo{x), the velocity field at a later time t is recovered from the 
potential ^{x,t) given by the relation 



exp 



dy 

exp 



1 , i^-y) 



2 \ 1 



*o(2/) 



2v V 2t 



(5) 



with $o(x) standing for the initial potential such that uo(x) = dx'^oix). The invicid Burgers equation 
corresponds to the limit v ^ Q. The solution is then given by solving a minimalization problem: 

u{x, t) = dx<^{x, t) with $(a;, t) = mm {^o{y) + ^ • (6) 

Outside shocks the minimum is reached for only one value of y, the solution ofthe equation uo{y^,)t = x—y^. 
The velocity is u{x,t) = = uo(j/^). It is effectively a local solution of the invicid Burgers equation since, 
by the minimum condition defining y^, we have u{x,t) = uq{x — tu{x,t)). A simple geometrical construction 
of the solution is described in refs.|0, ||]. For large t, y» coincides approximately with one of the local 
minima of ^o{y) and it practically does not change under small variations of x so that, in between the shocks, 
the velocity is approximately linear with the slope 

Shocks appear when the minimum is reached for two values yi and ?/2 of y. Let $1^2 — $0(2/1,2) be the 
value of the initial potential at these points. Then eq. ^ allows one to determine the velocity profil Us{x, t) 
around and inside the shocks at finite value of the viscosity ly by expressing exp [— -^^^(x, i)] as the sum of 
contributions from the two minima. One obtains: 

Us{x,t) = 7 (^2;- ^(2/1+2/2)^ - ^ tanh (^£^ (^a;-6<- ^(2/1 +2/2)^^ (7) 

where Ms = 2/1 ^ 2/2 > and = '^\Zy^ ■ In the invicid limit — )■ 0, this becomes 

us{x,t)l^, = 6 + ^^1^ - ^ (^9ix - x,M) e{x,{t) - x)) (8) 

where Xs{t) — + ^(2/1 + 2/2) is the time t position of the shock which moves with the velocity and 
follows a Lagrangian trajectory. 9{x) is the step function. The values of the velocity on the two sides of the 
shock are: 

uf = u,{xs±Q)=^sTYt ^^'^ 
so that ^ is the amplitude of the shock. 

The presence of shocks is at the origin of universal features which are independent of the details of the 
statistics. They may be analyzed by looking at fields localized on the shocks. By definition, these fields may 
be represented for any realization as: 

Og{x,t)^ J2 gi^s,^^s) S{x-x,{t)) (10) 

shocks 

where the sum is over the shocks with Xs(t) denoting the position of the shock, its velocity and ^ its 
amplitude. These fields are labeled by functions of and /is- By using the velocity profile (0) inside and 
around the shocks, we may map fields defined in terms of the velocity u{x,t) into the shock representation. 
For example, the shifted derivative of the velocity field {dxu{x,t) — j) is for large t localized on the shocks 
since away from shocks, u{x,t) — ^^—^ with ?/* almost independent of x. More generally, the generating 
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function (^dx — j) exp [\u{x, t)] also becomes localized on the shocks for large t. Using the velocity profiles 
(|^) or directly (|[), one finds its shock representation: 

dx - j) e^ = -2^e^«= sinh(^) Six - x,{t)) . (If) 



Note that this differs from the products X{dxu{x,t) — 1) e'^"'^^"'*) which are evaluated using the fact that 



t ' 

the velocity on the two sides of a shock are = ■Cs T — ' 



2t ■ 



A (dMx, t) - e^-(-±o-*) = - A^ ^ exp [a(C. T ^)] S{x - x^t)) . (f 2) 

Another example of a field localized on shocks is provided by the dissipation field e{x,t) defined by 

e{x,t) = lim i'(dxu)'^. As is well known, e{x), which is naively zero due to the prefactor v in its definition, 

1/— >o 

is actually a non-trivial field since (dxu)^ is singular in the invicid limit, fntegrating v(dxUs)^ around the 
shock one obtains in the limit ^ the contribution j^i^)^- Hence the shock representation of e(x) is: 

s 

More generally one finds, using again the velocity profile (|7|), the shock representation of the generating 
function e(a;, i) e'*' "*^^'*^. Namely: 

e,(x,t)^<x,Oe^"(-*) = 2A"3^e^«= cosh(^) - sinh(^)) - . (f3) 

Note that for v 0, the Burgers equation (|I|) implies that 

= (a* + xdxjdx-xty{d^u)'^c^-^ = (dt + xdxjdx + x^iy{dxuf-i,d^^e^". (u) 

Since S^e"^" has a (distributional) limit when — > 0, we may expect that at = 

(^dt + Xdxjdx)e^" + X-'ex = (f5) 



encoding the invicid version of the Burgers equation. Indeed, eq. ( |f S]) may be verified directly at = by 
computing (dt + Xdxjdx^ e^'' with the use of the limiting profile (^). 

Comparison of eq. (^3|) with eqs. ( p^JI^ ) yields an alternative representation of the dissipation field in 
terms of the velocity field: 

exix,t) = ^A-^ (a:,w(a;,<)) (2e^"(^^*) -e^"(^+°'*) -e^"(^-°'*)) (f6) 
^ iA-2 (e^ «(-'*)) (2a,M(x, t) - dxu{x + 0, t) - dxuix - 0, i)) 

Eq. (^6|) is an extension of the well-known formula e{x) — \inidi [u{x) — u(x + l)]^ . As expected and 

manifest in eq. (|I^) , the dissipation field is located on the discontinuity of the derivative of the velocity 
field. Eq. (p^ ) does not coincide with the operator product expansion suggested in ^ for the forced Burgers 
turbulence and expressing as a combination of e'^ " and i9a;e^". 

Fields localized on shocks form a closed algebra. When shocks are diluted, these operators are expected 
to satisfy a simple operator product expansion: 

Of{x,t)-Og(y,t) = S{x-y) Ofg{x,t) + regular. 

The contact term 5(x — y) in this operator product expansion arises from the coinciding shocks in the double 
sum representing the product operator. As an application, let us present an argument showing that the 
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structure functions scale linearly in x. Indeed, using the representation (^l]) for the operator (dx — A) e'^"(^'*) 
at i = 1, one finds: 

(9, - A)' (e^ =<5(x) (O^) + regular 

with Oip{x) = 2 ^j,(cosh(A/is) — 1) S{x — Xs). By integrating, this implies: 

^gAKx)-«(o))^ = 1 + a(A) X + \x\ + o{\x\) (17) 

with a(A) = — a(— A). Eq. ( p^ implies that, at short distances, (('«(a;) — u(0))") is proportional to |a;| for 
even positive n: 

{{u{x)~u{Q)r) = {0,,.) \x\ + o{\x\) 

and it is consistent with the behavior proportional to x for odd n > 1, as the one holding for the 3- 
point function. In words, the anomalous scalings of the structure functions in Burgers turbulence are a 
simple echo of the shocks. They are universal (at least when shocks are diluted): only the amplitudes 
are statistics dependent. Some of these representations and formal manipulations also apply to the forced 
Burgers equation. 



3 Self-similar solutions 

Self-similar behavior such as in eq. (^) will be true at any time, i.e. not only asymptotically, if the initial 
correlation functions scale. Indeed, from the explicit solution ^ it immediately follows that (for the Gaussian 
initial data) 

Goisx) ^ s'^'' Go{x) =^ s^-^ u{sx, s^-^t) u{x,t) . (18) 

h is the dimension of the initial potential. Such scaling behavior corresponds to a characteristic length 
l{t) ^ . For this length to grow with time we must have h < 2. Shocks are expected to be dense for 
1 < /i< 2 and diluted for /i< 1 |, I). 

Demanding a self-similar behavior for the correlation functions imposes constraints on the correlation 
functions of the velocity field: 

{2- h)tdt + J2{x,dx^ + {I - h)X,d,^)]{l[e'''-(^''^'^) = 0. (19) 

j k 

One may construct self-similar solutions with scaling dimensions h for > /i > — 1 by generalizing the 
representation of Kida's asymptotic solution ||]. By construction, the velocity Uh{x,t) has the following 
form: 

Uh{x,t)^dx^h{x,t) with ^h{x, t) ^ mm (^cbj + ^ (20) 

where ((f>j,yj)j^z is a Poisson point process with intensity fh{(f>)d(f>dy. Recall that this means that the 
probability to find a point of this process in an infinitesimal cell centered at ((/>, y) is //i(0, y) d(j)dy and that 
such elementary events are independent. To assure the translation invariance, fh will depend only on (p. 
For any given realization, the velocity field ( pO| ) has an exact sawtooth profile with slope \. In this Ansatz 
all shocks are created at time t = Q. The later time evolution is then governed by the shock collisions: the 
biggests eating the smallests. 



Let us first show that eq. (20) is preserved by the evolution specified by the invicid Burgers equation. At 
a time t' = t + t > t, the velocity field Uh{x, t') is given by: 

Uhix, t + T)^dx min f $/,(?/, t) + l^L—lI 
V \ 2t 
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Inserting the expression ( |20| ) for ^{y,t) and commuting the minimization over y and over j's, we get: 



mm 



2t 2t 



dx min 



J V 2(i + T 

Next we determine the intensity fh{(j)) d(f)dy such that the solution ( pO| ) is self-similar with scaling di- 
mension h, i.e. s^^^ Uh{sx, s'^^^t) ^ u{x,t). Let us spell out this condition for the potential ^h{x,t). By 
definition, 



mm , ^ , , 



where — s ^(t>j and yj = s ^j/j. Since Uh{x, t) = 9j;<i>/i(x, i), demanding self-similarity amounts to require 
that s~^^h{sx, s^^^t) = ^h{x, t) — Cs with Cg a constant. This equality will be true in law if the intensity is 
such that: fh{4>) d(j)dy = fh{4' + C's) d(j)dy. Up to an irrelevent translation of the solutions of this equation 



are: 



= ( ^ '^'^ for > ^.^j^ -l<h<Q, (21) 

lo for</.<0 ^ ^ 

/o((/)) = exp[const.0] with h = Q. (22) 



Thus we showed that the representation ( P0|) of the velocity field in terms of the Poisson process {(l)j,yj) 
with intensity (pl| ) is (i) self-similar and (ii) preserved by the evolution. In particular the relations such as 
eqs. (|r^ ) will be satisfied. Moreover the invicid form (^) of the Burgers equation for each realization implies 
the Hopf equations for the correlators: 

5t + ^A,a,^ ^a,^.]([]e^'="(--*)) + (e,^ [] e^'="(--*)) = (23) 

j ^ k j k^j 

The time derivative may be eliminated from both equations leading to the fixed-time version of the Hopf 
equations. The case h = corresponds to Kida's asymptotic solution. As pointed out in |^ it generalizes 
the Gumbel class of extreme statistics. The case h ^ should, correspondingly, generalize the WeibuU class 
of extreme statistics. 

Let us describe in more details the case h = It corresponds to initial potential correlation functions 
homogeneous of degree —1, e.g. like the Dirac delta function. We denote by ^*{x,t) and u*{x,t) the corre- 
sponding potential and velocity field. The statistics of the Poisson point process {(f)j,yj) in the representation 
( pO| ) is specified by the intensity 

f*(^\ - / (l>d(j)dy for (/) > 0, 

^'^^ " \ for (/.<0. ^ ' 

D is a constant with dimension (length)^ x {time)^'^. Recall that the velocity field u*{x,t) is such that 
u*{x,t) = t^^/^u*{xt^'^/^ , 1). In other words, the characteristic length at time t is l*{t) — D^/^ t^l^ . 

The one-point function of the potential scales as t^'^l'^ and thus diverges at t = (it does not contribute 
to the one point function of the velocity which vanishes). The two-point function G*(x, t) of ^*{x, t) satisfies 
G*{x,t) = t~^/^G*(a;t~^/^, 1). Since, as we shall see, G*{x, 1) is smooth, regular at the origin and decreasing 
exponentially at infinity, at zero time G{x,t) becomes proportional to the Dirac delta function, 

lim G*(x, t) — const. S(x) . 

Once the velocity field has been parametrized in terms of Poisson processes as in eq. (^0|) , it is easy to 
compute any correlation functions. For example, the one-point generating function is 

(e^"*("'*^) = f d(t,dyVx{<l>,y) e^(--^')/* 
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with Vx{4>,v) denoting the probabiHty (density) of {4>,y) = {<l>j' lUj*) for j* minimizing ^'/'j + 2t^ ^ ) • 
To compute this probabihty, imagine dividing the ((/), ?/)-plane into elementary cehs of size d(j)dy with the 
probabiHty for a point {(j)j,yj) to be in the cell centered at (i^i', y) equal to f {(j))d(j)dy . Thus, the probabihty 
Vx{(f>, y) is equal to the product of the probability for the elementary cell centered around ((/>, y) to be occupied 
by a point of the process times the probability for the other cells around (0', y') s.t. <j)' + 21 < 0+ 2t 
to be empty (if this condition is violated, the cell around (0', y') may be either occupied or empty). Hence 

'Px{4>,y)d(l)dy = f{(t>)d(t>dy J| [l - x{4>\v' \4>,v)f{4>')d(t>' dy' 

d(j)' dy' 



where xi^t'' tV'] 4'Ty) is the characteristic function of the constraint 



< 



Taking the 



2t V I 2t 

continuum limit of infinitesimal cells and approximating the product over the cells around {<j)' ^y') by the 
exponential of a sum leads to the result: 



^Au*(2;,l)\ 



d<pdy e-^y /(0-i^)exp 



dz 



2 

w 



/ dX X{X coshX - sinhX) cxp 
Jo 



(25) 



We have set t = 1 and D — 1. The dependence on these parameters is restored by replacing A by XD^^^t~^/^ . 
One may check directly that the above expression for the 1-point function satisfies the identities (^9|) and 
(]2^). (c**" is regular around A = and its behavior at infinity is: 



const. lAj-i^/^ exp[-const. lAI^/'^] for A -> 00 . 



(26) 



This has to be compared with Kida's statistics for which (c*"' (^'i)) = exp[— const. A^]. The two-point 
functions can be computed similarly. For x > 0, 



^\{u'{x}-u*{-x))\ 



„2\x 



d(j)dy 



<t>>0 



(j) + 2xJa(0; y\ a;) Ja(0; -y; x) 



-Io{<l>;y;x)-Iol4>;-y;x) 



(27) 



with 



dz{(t) + -{x + yf 



\ J^zic^+^-ix + yf-lz^ 



where the domain of integration in both cases is V ~ { z \ z < y + x; z^ < 2(f> + {y + x)"^}. Eq. (|2^ ) may be 
used to check that for t 7^ the two-point function G{x,t) is smooth, fast decreasing at infinity and regular 
at the origin. There arc no difficulties, but not much motivations, to compute in the same way the higher 
point correlation functions. 



4 Comments and speculations 

We have constructed self-similar solutions of the decaying Burgers turbulence. It is natural to wonder if 
such statistics effectively describe the long time behavior of systems with smooth random initial data. One 
may construct such examples in a tautological way by taking as initial data the potential obtained 
from the self-similar Ansatz ( pC| ) at small but non zero time: ^q{x) = ^hix, rj) with rj ^ 0. By construction, 
it defines a smooth initial statistics which will have a large time asymptotics given by ^h{x,t). Of course, 
this initial statistics is not Gaussian. This shows, however, that the basin of attraction of the self-similar 
solution ( po| ) is not totally empty. 

In the case h — the two-point function of ^q{x) tends to the Dirac delta function 5{x — y) when 77 — > 0. 
So the initial statistics ^q{x) may be thought of as a way to regularize an initial potential with d{x — y) 
two-point correlation function. Note that if we replace 6{x — y) by a smooth cut-off dependent function 
A(x — y), the values at the origin A(0) diverge with the cut-off and Kida's asymptotic regime disappears in 
the limit since its characteristics length diverges. 
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The question is then whether there exist analogues of eq. (Q) but with different respective scaling between 
X and t corresponding to non-zero values of h. For example one may inquire about the existence of the limit 

Naively, it corresponds to a limiting initial potential with S{x ~ y) correlation function. Indeed, upon 
assuming that min„ and lim commute, it follows from eq. (|) that: 



lim -L $ f ^, = lim -L min ( $(y) + £5/2(^1^^! 

e^O el/2 I g ' g5/2 / el/2 y \ 2t 



mm[n*iy) + K^] (29) 



y 



2t 



where 'i>Q*(^) = 1™ 7T72 ^ (f ) the rescaled initial potential with the two point function G**{x) = 

lim i-G (1) — D d{x) with D = J dxG{x). It has dimension h — — i. Such a limit ( |2^ ) would correspond to 

an intermediate regime with characteristics length l*{t) ~ t^/^ smaller than Kida's length Ixit) ^ t^^^. The 
problem, however, is that the realizations of the white noise $q* are distributional and the last expression in 
( |29| ) is ill-defined. Clearly a finer analysis is required to decipher cases in which a limit of the type ( p8| ) leads 
to a self-similar asymptotic distribution with h — — i as the one constructed above. In fact, we expect the 
latter to appear if the initial potentials have values uniformly bounded below, in analogy with the extreme 
statistics problem. 

Identical constructions, arguments and speculations could be done in higher dimensions. For example, 
in dimension d the delta function has dimension /i = — ^, and there exists a self-similar solution with this 
scaling dimension. It corresponds to the characteristic length l{t) ~ t" with a — 

Finally, we feel that it could be worth-while to adapt the renormalization group techniques to analyze 
this type of large time behavior. 
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